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O ■ Comparing the gravitational acceleration induced on the Local Group by optical (SSRS2; 
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da Costa et al 1994, 1998) and IRAS (1.2 Jy; Fisher et al 1995 and 0.6 Jy - QDOT; Rowan- 
Robinson et al 1990) galaxies we estimate, within the framework of linear theory, their 
o relative bias factor. Using both IRAS samples we find boi ~ 1.14 — 1.2, slightly lower than 

O \ in Willmer, Da Costa & Pellegrini (1998) who use a £(r) approach. 

CO 



Using linear perturbation theory one can relate the gravitational acceleration induced on the 
observer by the surrounding mass distribution to her/his peculiar velocity: 



u(r) = irz^ J 6 ^w dr = ir DW 



where 



1 N 1 

4vrn ^ farf 



with 0(r) = (ra)" 1 J i ma ^ T .^ $(L) dL, where 3>(L) is the luminosity function of the objects 
under study and L min (r) = 47rr 2 Si im , with S\ im being the flux limit of the sample under 
study. 

Different classes of extragalactic objects (QSOs, AGNs, galaxies, clusters of galaxies) trace 
differently the underlying matter distribution, usually assumed to be linearly related via 
(^P/p)traccr = & traccr,m (<Wp) mass ( cf - Kaiser 1984). Although this complicates the use of 
dipole as an estimator of the cosmo logical density parameter, Q , it does allow one in principle 
to study the relative bias displayed by such objects (Kolokotronis et al 1996; Plionis 1995). 
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Using different tracers, ie. SSRS2 optical and IRAS galaxies, to determine the Local Group 
acceleration (dipole) we can write: u(r) = Q®' 6 D (r)/b = f^ ,6 Di(r)/&j and therefore we 
can obtain an estimate of their relative bias factor from: 



boi = 



Do(r) 
Dr(r) 



r(r). 



A statistically more reliable approach is to use the differential dipole, estimated in equal 
volume shells, to fit boi according to: 

2 _ ^ s (D 0) i — 6ojD/,i — Ci) 2 

where C is the zero-point offset of the relation and a are the corresponding shot-noise errors 
estimated according to Strauss et al. (1992). 



6000 



J I I I I I I L 



6000 




50 100 150 

R (h" 1 Mpc) 



_l I I I I I I I L 



_l_ 



50 100 150 

R (h" 1 Mpc) 





1.8 




1.6 


s- 




o 






1.4 


o 


CO 




Si-i 




t/I 


1.2 


CO 




2 






1 




0.8 



i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 
Ogopt/g IRAS 

differential Dipole fit J 



tf 




_L 



b 0I =1.16 

_l I I I I I L 



50 100 

R (h" 1 Mpc) 



150 



1.8 
1.6 h 
1.4 



n 1.2 

CO 



0.8 - 



T 



T — I — I — I — I — I - 




b 0I =1.18 
L 



J I I I I I I I I I I I I L 



50 100 150 

R (tr 1 Mpc) 



2 



In the figure and in table 1 we present the results of both methods used to estimate the 
optical (SSRS2) to infrared (IRAS 1.2 Jy and QDOT) bias factors. 

Table 1: Optical to Infrared galaxy bias factors from differential dipole fit. 
samples boi C \ 2 d.f. 

SSRS2-QDOT 1.16 ±0.13 -80 ±70 8.7 6 

SSRS2-IRAS 1.2 Jy 1.18 ±0.06 -30 ±30 15.7 6 
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We construct the X-ray luminosity (LF) for different classes of galaxies (Seyferts, el- 
lipticals, star- forming galaxies and Liners), by convolving the optical LF of the Ho et al. 
spectroscopic sample of nearby galaxies with the corresponding L x /L relations from the 
Fabbiano et al. X-ray atlas of galaxies. From the derived LF we can easily assess the con- 
tribution of galaxies to the X-ray background. The Seyferts and Liners make the largest 
contribution (~ 40%) assuming no evolution while the contribution of star-forming galaxies 
is much smaller (~ 5%). 

The Sample & Method 

We have used the spectroscopic sample of nearby galaxies (B<12.5) of Ho et al. (1995). The 
great advantage of this sample is that there is excellent spectroscopic infomation (high signal- 
to-noise, medium resolution) available and thus bona-fide spectroscopic identifications exist 
for all (~ 500) galaxies in the sample. Hence, we can construct the X-ray LF separately for 
Seyferts, star-forming galaxies (HII), ellipticals and Liners instead of simply dividing them 
to early-type and late-type according to their morphology. The majority of galaxies in the 
Ho et al. sample are HII (50%), Liners are (30%), Seyferts (13%) while 15% of galaxies 
present no emission lines and thus can be classified as 'early-type' galaxies or 'ellipticals'on 
the basis of their spectra rather than their morphology (Ho et al. 1997). We first derive 
the optical LF for different classes of objects (Seyferts, HII, Liners and no-emission-line or 
early-type galaxies) using the classical 1/V max method. Next, we derive the L x /Lb relation 
for the different subclasses using the EINSTEIN X-ray fluxes (0.5-4.0 keV) from the Fabbiano 
et al. (1992) X-ray atlas of galaxies. There are 164 entries (95 detections and 69 3a upper 
limits) of Ho et al. galaxies in the Fabbiano atlas. Finally, in order to derive the X-ray LF 
we convolve the optical LF with the L x /Lb relation: 
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$(L X ) = J $(L B )(/)(L x \L B )dL B 

(eg Avni & Tananbaum 1986), where <j)(L x \L B ) is the conditional probability function and 
can be expressed as a Gaussian around the mean L x value for a given L B (we use H Q = 100 
throughout). 1 1 1 1 1 1 1 1 



44 



42 



-1 



36 



•O • CT^i 



I i i i I i i i I 

42 44 46 

lo 



Although, our method provides only an indirect way of deriving the X-ray LF, it is currently 
the only feasible method at least for some classes of galaxies. Unfortunately, the galaxies in 
X-rays are faint (apart from Seyferts) and thus we cannot yet obtain large X-ray selected 
galaxy samples neither in deep X-ray surveys nor in the ROSAT all-sky survey. For example 
we note that there are only 5 HII galaxies from the Ho et al. sample detected by the RASS 
(Zezas et al. 1998). 

The Results 

In the figure we plot the X-ray luminosities vs. the Blue luminosities (log 10 (L x ) vs log 10 (L s )), 
for different subclasses (open symbols represent upper limits while filled symbols represent 
detections; the color coding is as in the table). 

The X-ray LF is plotted on Fig. 2. It is clear that although the HII galaxies are the most 
numerous, due to their low they contribute much less X-ray luminosity per h -3 Mpc 3 
compared to Seyfert galaxies. Indeed, the HII emissivity is J ~ 1 x 10 38 erg s _1 h 3 Mpc -3 
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Table 2: The best-fit $ x results for a Schechter optical luminosity function 
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while the Seyfert is J ~ 8 x 10 38 erg s _1 h 3 Mpc -3 . The total galaxy emissivity is comparable 
to the emissivity of local galaxies in the 2-10 keV band derived by Lahav et al. (1993) from 
the cross-correlation of optical galaxy catalogues with the fluctuations of the hard X-ray 
background. 
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We estimate the dipole of the diffuse 1.5 keV X-ray background from the ROSAT all-sky 
survey map of Snowden et al (1995). We first subtract the diffuse Galactic emission by fitting 
to the data a finite disk model, following Iwan et al (1982). We further exclude regions of low 
galactic latitudes, of local X-ray emission (eg the North Polar Spur) and model them using 
two different methods. We find that the ROSAT X-ray background dipole points towards 
(1, b) « (288°, 25°) ± 19° within ~ 30° of the CMB; its direction is also in good agreement 
with the HEAO-1 X-ray dipole at harder energies. The normalised amplitude of the ROSAT 
XRB dipole is ~ 1.6%. 



Subtracting the Diffuse Galactic Emission 

A major problem in extracting the X-ray background dipole at soft X-rays is the large 
contribution of our Galaxy at such energies (cf. Kneissl et al 1998). We attempt to model the 
diffuse Galactic component using a finite radius disk with an exponential scale height, (Iwan 
et al. 1982) which provides a good description of the Galactic component at harder energies 
(2-60 keV). The fit to the model is performed by excluding the regions of the most apparent 
extended Galactic emission features: the North Polar Spur and the |6| < 20° strip as well 
as most apparent "local" extragalactic features; a 4° radius region around the Virgo cluster 
(/, b « 287°, 75°) and a 10° radius region around the Magellanic clouds (I, b « 278°, -32°). 

We obtain a fraction of the total X-ray emission which is due to the Galaxy consistent 
with ~ 25% while the disk scale height and disk radius where found to be 16 & 27 kpc 
respectively. 
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Dipole Estimation 



After excluding from the ROSAT counts our best fit Galactic model and after masking, using 
either of two methods; homogeneous filling procedure or a spherical harmonic extrapolation 
procedure, the following regions: (a) The Galactic plane, with |6| < 20° or 30°, (b) the 
area dominated by the Galactic bulge and the North Polar Spur (ie., —40° < b < 75° and 
300° < / < 30°) , (c) the Large Maggelanic Clouds (ie., an area of 10° radius centred on 
I, b ~ 278°, —32°), we measure the dipole by weighing the unit directional vector pointing to 
each 40 2 arcmin 2 ROSAT cell with the X-ray intensity Cj of that cell. We normalize the dipole 
by the monopole term (the mean X-ray intensity over the sky): T> = \~D\/M = J2Ci*i/ Z)Ct 

Dipole Results 

When using the raw ROSAT data, the dipole points towards the Galactic centre (in agree- 
ment with the analysis of Kneissl et al 1997). However, when we exclude both the Galaxy 
and the North Polar Spur, the measured dipole is in much better directional agreement with 
the CMB dipole. For the homogeneous filling method we find S8 cm b ~ 20° while for the 
spherical harmonic method the misalignment angle is larger, 56 cm b ~ 51°. 
However, it should be expected that many Galactic sources, probably dominating the higher 
ROSAT counts, are still present in the data and could affect the behaviour of the estimated 
extragalactic XRB dipole. Excluding the highest ~ 3% of the counts (C up ^ 140) we find that 
both methods used to model the masked areas agree and the XRB-CMB dipole misalignment 
angle is reduced significantly 

The interpretation that the high intensity cells are associated with Galactic sources is sup- 
ported by the fact that when we include these few cells the resulting dipole direction moves 
towards the Galactic centre. 

Therefore, taking into account the variations of our results due to (a) the uncertainties of 
the Galactic model subtracted from the raw counts, (b) to the different methods used to 
mask the excluded sky regions and (c) to the different galactic latitude limits, we conclude 
that the extragalactic ROSAT dipole has: 

2>xrb,rosat « 0.016 ± 0.008 (I, b) « (288°, 25°) ± 19° 

Our results are consistent with the HEAO-1 (2-10 keV) dipole (Shafer & Fabian 1983) which 
points in a similar direction (282°, 30°) but has a lower amplitude: £>xrb,heao-i ~ 0.005. 
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We investigate the clustering properties of a new homogeneous sample of 2000 extra- 
galactic sources selected from the ROSAT XRT/PSPC All-Sky survey (Voges et al 1996). 
Assuming that this sample is dominated by QSO's we find that the characteristic depth of 
the sample is ~ 410 h~ l Mpc. We find a significant angular correlation function, between 
~ 2° and 15° which roughly corresponds to a spatial correlation length of r Q ~ 15.5 ±4.5 h~ x 
Mpc, roughly consistent with that of optically selected QSO's at (z) ~ 1.5. 

Method 

The original RASS sample contains 18811 sources which include stars and extragalactic 
sources over the entire 0.1-2 keV energy range. To exclude the stellar objects we cross 
correlated this sample with star catalogs (SAO, GSC, RASSOB, RASSWD, XRBCAT and 
CVCAT). Furthermore, we have cross-correlated the remaining sources (for b > 0°) with 
the Hamburg identifications (Bade et al 1998) to find that most of them, after excluding 
those with extension flag > 40, are identified as QSO's. In order to produce a homogeneous 
sample in exposure time we have omitted sources with count rate < 0.1 which results in a 
homogeneous sample over 92% of the sky. Finally, we exclude regions heavily affected by 
Galactic absorption, ie., |fe| < 30° as well as 5 < —30° (see below). Our final catalogue 
contains 2000 extragalactic sources. 

We further need to estimate the various selection functions that could affect the clustering 
properties of our survey. For example the significant absorption due to the diffuse Galactic 
Neutral Hydrogen could artificially enhance the 2-point correlation function and also intro- 
duce large-scale modulations. Therefore, we have cross correlated our RASS catalogue with 
the recent Leiden/Dwingeloo Atlas of Galactic Neutral Hydrogen (Hartmann & Burton 1997) 
which covers regions with 5 > —30°, to derive the RASS surface density as a function of Nh, 
which is used to produce a random catalogue with similar absorption selection function. 
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Results 



The correlation function estimator we use is: w{9) = fNdd/Ndr — 1, where the normalization 
factor is / = 2 x N r /N . Assuming a 2-p correlation function of the form w{6) = {6/6 ) h we 
find b = -0.9 ± 0.15 and O = 0.08° ± 0.05°. 

A preliminary calculation using the average flux limit (~ 10 -12 erg/sec cm 2 ) of the RASS 
over the area used and a QSO characteristic luminosity of ~ 2 x 10 (Boyle et al 1993) 
gives a characteristic depth of D* ~ 410 h" 1 Mpc. Using Limber's equation to relate the 
angular and spatial 2-point correlation function with 7 ~ 2 we find a correlation length of 
r D ~ 15.5±4.5 h~ l Mpc, which is compatible only with comoving QSO clustering (cf. Shanks 
k Boyle 1994 and La Franca et al 1998 which find r ~ 10 ± 2.5 h" 1 Mpc at (z) ~ 1.5 — 1.8) 




0.001 
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We study the performance of two different methods used to define cluster shapes with 
final aim to study the projected and intrinsic shapes of APM clusters (Dalton et al 1997). 
The I s * method defines the cluster ellipticity by fitting ellipses to the individual galaxy 
distribution as a function of radius from the cluster centre. The 2 nd method is based on 
smoothing the discrete galaxy distribution and then fits an ellipse using all cells that fall 
above an overdensity threshold. This latter method is free of the known aperture bias which 
tends to artificially sphericalize clusters. Using Monte-Carlo simulations we have studied 
the performance of both methods in the presence of the expected galaxy background at the 
different distances traced by the APM clusters. 

The Method 

In order to estimate the projected cluster shape we diagonalize the inertia tensor (det(L_j - 
A 2 M 2 ) = 0) where M 2 is the 2x2 unit matrix. The eigenvalues (a, b) with (a > b) define the 
ellipticity of the configuration under study: e = 1 — b/a (cf. Plionis, Barrow & Frenk 1991). 

Initially all galaxy positions are transformed to the coordinate system of each cluster by 
: Xi = (RA g — RAd) x cos(Sd) and = 5 g — 5 C [. Then we use the following procedures: 
I s * Method: All galaxies within an initial small radius are used to define the initial value of 
the cluster shape parameters. Then the next nearest galaxy is added to the initial group and 
the shape is recalculated. This method, although straight forward, suffers from the fact that 
implicitly we are assuming a spherical aperture within which the cluster shape parameters 
are estimated (cf. Binggeli 1982). 

2 nd Method: The discrete galaxy distribution is smoothed using a Gaussian kernel. Then 
all cells that have a density above some threshold are used to define the moments of inertia 
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with weight w { = p. t - (p(z)) / (p(z)) where (p(z)) = ImZu*) ®( M ) dM with $ ( M ) the APM 
luminosity function (with z evolution) from Maddox, Efstathiou & Sutherland (1996) and 
M m i n {z) = mu m — 5 log(r) — 25 — 3z. This method is free of the aperture bias and we found 
that it performs significantly better than the previous method. 



A major problem in reliably determining cluster shapes is the significant galaxy background 
which contaminates the cluster galaxy counts, especially in deep galaxy catalogues like the 
APM which has mn m = 20.5. In order to assess the effects of such contamination on the 
determination of cluster shapes we have performed large sets of Monte-Carlo simulations in 
which we compared the determination of cluster shapes, defined to have the King's profile 
with core radii < 0.1 h" 1 Mpc (cf. Girardi et al 1998), with and without the expected 
background contamination. The background at each cluster distance was estimated by: 



where 9 l is the ancular radius of the cluster. In the figure we show an example of a simulated 
cluster at 200 h~ x Mpc with e = 0.5. The performance of the two methods, described above, 
is shown in the same colour coding as in the text. It is evident that the 2 nd method performs 
equally well in the presence or not of the galaxy background. However, the presence of 
substructure can affect unexpectedly the performance of this method as well. 

We are at the process of quantifying biases related to the above methods and applying 
our shape determination procedure to the real APM data. 
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SIMULATED CLUSTER at 200 h" 1 Mpc with e = 0.5 & King's profile 
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